ARL-TR-8009  •  Apr  2017 


ARL 

US  Army  Research  Laboratory 


Tunable  Solid-State  Quantum  Memory  Using 
Rare-Earth-Ion-Doped  Crystal,  Nd3+:GaN 


by  S  Rudin,  BC  Connelly,  VS  Malinovsky,  RW  Enck,  and 
GD  Metcalfe 


Approved  for  public  release;  distribution  unlimited. 


NOTICES 

Disclaimers 

The  findings  in  this  report  are  not  to  be  construed  as  an  official  Department  of  the 
Army  position  unless  so  designated  by  other  authorized  documents. 

Citation  of  manufacturer’s  or  trade  names  does  not  constitute  an  official 
endorsement  or  approval  of  the  use  thereof. 

Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return  it  to  the  originator. 


ARL-TR-8009  •  Apr  2017 


ARL 

US  Army  Research  Laboratory 


Tunable  Solid-State  Quantum  Memory  Using 
Rare-Earth-Ion-Doped  Crystal,  Nd3+:GaN 


by  S  Rudin,  BC  Connelly,  VS  Malinovsky,  RW  Enck,  and 
GD  Metcalfe 

Sensors  and  Electron  Devices  Directorate,  ARL 


Approved  for  public  release;  distribution  unlimited. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  the  collection  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  the 
burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302. 
Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently 
valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  To) 

April  2017  DRI  Report  October  201 4-October  2016 

4.  TITLE  AND  SUBTITLE  5a.  CONTRACT  NUMBER 

Tunable  Solid-State  Quantum  Memory  Using  Rare-Earth-Ion-Doped  Crystal, 

Nd3+:GaN  5b.  grant  number 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S)  5d.  PROJECT  NUMBER 

S  Rudin,  BC  Connelly,  VS  Malinovsky,  RW  Enck,  and  GD  Metcalfe  DRI14-SE-002 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PERFORMING  ORGANIZATION  REPORT  NUMBER 

US  Army  Research  Laboratory 

ATTN:  RDRL-SEE-I  ARL-TR-8009 

2800  Powder  Mill  Road 
Adelphi,  MD  20783-1138 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES)  10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT  NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

The  research  objective  of  this  Director’s  Research  Initiative  was  to  work  on  developing  solid-state  quantum  memory  using  cryogenically 
cooled  rare-earth-ion-doped  crystal,  Nd3+:GaN.  The  samples  were  grown  using  the  US  Army  Research  Laboratory’s  molecular  beam  epitaxy 
(MBE)  capabilities.  The  approach  is  based  on  the  coherent  transient  interaction  between  light  and  the  ions  in  a  semiconductor.  We  investigated 
the  energy  level  structure  of  the  neodymium  (Nd)  ions  embedded  into  the  semiconductor  crystal  lattice,  gallium  nitride  (GaN),  and  developed 
a  quantum  description  of  the  2-  and  3-pulse  photon  echo  effect.  Epitaxial  GaN:Nd  was  grown  by  plasma-assisted  molecular  beam  epitaxy  in 
a  modular  Gen  II  reactor  using  liquid  gallium,  solid  Nd,  and  a  nitrogen  plasma.  The  photoluminescence  (PL)  spectrum  from  the  lowest  4F3/2 
multiplet  to  all  4U/2  multiplets  was  measured  by  a  free- space  variable  amplifier  silicon  diode.  A  demonstration  of  the  coherence  properties 
requires  a  high-quality  sample  with  a  narrow  homogeneous  broadening.  For  this  reason,  the  PL  linewidth  was  characterized  for  various  growth 
conditions.  The  experiments  were  carried  out  to  observe  coherent  effects,  and  various  experimental  setups  are  described.  In  parallel  with 
experiments  on  the  coherence  properties  of  Nd3+:GaN,  we  also  investigated  the  effect  of  a  varying  electric  field  on  the  Nd  PL  spectrum  to 
provide  a  tunable  memory.  To  vary  the  applied  field,  we  designed  and  grew  a  series  of  Nd-doped  GaN  p-i-n  structures,  strain-balanced 
superlattice  structures,  and  periodic  unbalanced  structures. 


15.  SUBJECT  TERMS 

solid-state  quantum  memory,  neodymium,  2-level  system,  Stark  shift,  hetero structure 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF 

ABSTRACT 

UU 

18.  NUMBER 

OF 

PAGES 

44 

19a.  NAME  OF  RESPONSIBLE  PERSON 

S  Rudin 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

19b.  TELEPHONE  NUMBER  (Include  area  code) 

Unclassified 

Unclassified 

Unclassified 

301-394-0206 

Standard  Form  298  (Rev.  8/98) 
Prescribed  by  ANSI  Std.  Z39.18 


11 


Contents 


List  of  Figures  v 

1.  Summary  1 

2.  Introduction  1 

3.  Approach  to  Achieving  Research  Objectives  2 

3.1  Quantum  Memory  Based  on  2-Pulse  Photon  Echo  3 

3.2  Quantum  Memory  Based  on  Controlled  Reversible  Inhomogeneous 

Broadening  (CRIB)  4 

4.  Electronic  States  in  Rare  Earth  Ions  6 

5.  Review  of  Previous  Growth  of  Neodymium  (Nd)-Doped  Gallium 

Nitride  (GaN)  7 

6.  Theoretical  Analyses  of  2-  and  3-Level  Models  and  Design 

Parameters  8 

6.1  2- Level  System  8 

6.2  Bloch  Equations  11 

6.3  3-Level  System  and  Fast  Nonradiative  Decay  of  the  Top  Level  11 

6.4  3-Level  System  and  Fast  Nonradiative  Decay  of  the  Intermediate 

Level  13 

7.  Steady-State  Population  14 

8.  Pulsed  Excitation  of  the  2-Level  System  17 

9.  Experimental  Work  for  Nd3+:GaN  Director's  Research  Initiative 

(DRI)  20 

9.1  Growth  of  Nd3+:GaN  20 

9.2  Continuous-Wave  (CW)  Steady-State  Saturation  Photoluminescence 

(PL)  Setup  and  Experiment  21 

Approved  for  public  release;  distribution  unlimited. 

iii 


9.3  Characterization  of  Transition  Linewidth  23 

9.4  Time-Gated  Photoluminescence  (PL)  Setup  and  Experiment  for  Rabi 

Oscillations  and  Photon  Echo  Measurements  24 

9.5  Characterization  of  Field-Shifting  of  Transition  Wavelength  27 

10.  Conclusion  31 

11.  References  32 

List  of  Symbols,  Acronyms,  and  Abbreviations  34 

Distribution  List  36 


Approved  for  public  release;  distribution  unlimited. 


List  of  Figures 


Fig.  1  Evolution  of  the  Bloch  vector  during  the  2-pulse  photon-echo  process, 
a)  tc/2  pulse  rotates  the  Bloch  vector  from  the  negative  w  axis  along 
the  u  direction,  b)  The  individual  Bloch  vectors  precess  in  the  uv  plane 
and  dephase  due  to  the  inhomogeneous  broadening  of  the  transition,  c) 
At  time  t  =  x,  a  re-pulse  rotates  all  Bloch  vectors  around  the  v  axis,  d) 
The  Bloch  vectors  realign  and  build  up  a  macroscopic  coherence  that 


results  in  a  photon-echo  emission  at  time  t  =  2x . 3 

Fig.  2  a)  Stimulated  photon-echo  schematic,  b)  Generalization  of  the  3-pulse 
echo  to  data  storage . 5 

Fig.  3  Quantum  state  storage  based  on  CRIB  in  crystals  with  an  ensemble  of 
ions  having  a  broad,  inhomogeneously  broadened  absorption  line . 5 

Fig.  4  Low-temperature  (~30  K)  PL  from  900  to  1500  nm  for  GaN:Nd 
pumped  at  836  nm  showing  the  characteristic  Nd  transitions. 
(Reprinted  from  Readinger  et  al.13) . 8 

Fig.  5  Schematic  of  a  2-level  system . 9 


Fig.  6  Schematic  of  a  3-level  system  with  nonradiative  decay  of  the  top  level 


. 11 

Fig.  7  Schematic  of  a  3-level  system  with  nonradiative  decay  of  the 

intermediate  level . 13 

Fig.  8  Steady-state  population  in  2-level  system,  Eq.  8 . 15 

Fig.  9  Useful  dependencies  to  analyze  the  steady-state  population  as  a 

function  of  the  field  amplitude . 15 

Fig.  10  Steady-state  population,  see  Eqs.  8  and  15 . 16 

Fig.  1 1  Useful  dependencies  to  analyze  the  steady-state  population  as  a 

function  of  the  field  amplitude . 16 

Fig.  12  Steady-state  population,  see  Eqs.  8  and  20 . 17 

Fig.  13  Useful  dependencies  to  analyze  the  steady-state  population  as  a 

function  of  the  field  amplitude . 17 

Fig.  14  Integrated  PL  and  the  excited-state  population  vs.  time  for  the  limiting 
case  of  constant  Rabi  frequency . 19 

Fig.  15  Integrated  PL  and  the  excited-state  population  vs.  time  for  the  pulse 

excitation . 20 

Fig.  16  Rabi  frequency  vs.  time  for  the  pulse  excitation . 20 

Fig.  17  Simplified  energy  diagram  of  the  manifolds  in  the  Nd3+  ion  in  GaN. 

(Reproduced  from:  Metcalfe  G,  Redinger  E,  Shen  H,  Woodward  N, 


Dieroff  V,  Wraback  M.  Crystal-field  split  levels  of  Nd3+  ions  in  GaN 
measured  by  luminescence  spectroscopy.  J  Appl  Phys. 
2009;105:053101.,  with  the  permission  of  AIP  Publishing) . 21 


Approved  for  public  release;  distribution  unlimited. 

V 


Fig.  18  Left  panel:  Experimental  setup  for  CW  PL  and  time-gated  PL.  Right 
panel:  Emission  spectrum  from  the  lower  4F3/2  state  to  the  4L/2  ground 
multiplet;  the  arrow  indicates  the  transition  between  the  ground  levels 
of  each  state  at  916.2  nm,  which  is  the  strongest . 22 

Fig.  19  Left  panel:  Saturation  of  the  916.2-nm  emission  line  using  836-nm 
excitation  with  increased  incident  laser  electric  field  strength.  Right 
panel:  First  and  second  derivative  of  the  best  fit  to  the  saturation  data 
in  the  left  panel.  The  value  at  Y  =  0  of  the  second  derivative  is  related 


to  the  dipole  moment  of  the  transition . 23 

Fig.  20  PL  of  Nd3+:GaN  samples  with  varying  doping  levels.  Right  panel: 

Zoomed-in  plot  of  transition  near  916  nm . 24 


Fig.  21  Left  Panel:  Timing  diagram  for  the  measurement  of  Rabi  oscillations. 

The  gating  pulse  removed  the  in-pulse  PL  signal  and  laser  scattering. 
Middle  Panel:  Sample  time-resolved  PL  emission  signals  for  a 
sampling  of  pulse  widths.  Right  Panel:  Integrated  PL  signal  as  a 
function  of  gating  pulse  width  shows  only  a  linear  increase  in  signal. 25 

Fig.  22  Measurement  of  T\  =  150  ps  in  Nd3+:GaN . 27 

Fig.  23  Electric  field  vs.  GaN  well/  AlxGa(i-X)N  barrier  thickness  for  strain- 

balanced  superlattice  (SBSL)  structures  with  GaN  wells/  Alo.8Gao.2N 
barriers  (orange  line)  and  periodic  unbalanced  (PUB)  structures  with 
GaN  wells  and  Al0.2Ga0.sN  barriers  (green  line) . 28 

Fig.  24  SIMS  measurement  on  the  60-nm-thick  well/barrier  SBSL  sample...  28 

Fig.  25  Nd3+:GaN  PL  spectra  of  an  approximately  zero-field  bulk  sample  and 
high-field  SBSL  samples  (expected  electric  field  strength  given  in 
legend) . 29 

Fig.  26  TEM  images  of  sample  1521,  the  40-nm-thick  well/barrier  SBSL 

sample,  with  a  100-nm  scale  bar  (left  panel)  and  a  20-nm  scale  bar 
(right  panel) . 29 

Fig.  27  Electron  holography  phase  image  (left  panel)  with  area  for  line  profile 
(right  panel)  outlined.  Brighter  areas  are  GaN,  and  darker  areas  are 
AlGaN . 30 


Approved  for  public  release;  distribution  unlimited. 


1.  Summary 


The  research  objective  of  this  Director’s  Research  Initiative  was  to  work  on 
developing  solid-state  quantum  memory  using  cryogenically  cooled  rare-earth-ion- 
doped  crystal,  Nd3+:GaN.  The  samples  were  grown  using  the  US  Army  Research 
Laboratory’s  molecular  beam  epitaxy  (MBE)  capabilities.  The  approach  is  based 
on  the  coherent  transient  interaction  between  light  and  the  ions  in  a  semiconductor. 
We  investigated  the  energy  level  structure  of  the  neodymium  (Nd)  ions  embedded 
into  the  semiconductor  crystal  lattice,  gallium  nitride  (GaN),  and  developed  a 
quantum  description  of  the  2-  and  3-pulse  photon  echo  effect.  Epitaxial  GaN:Nd 
was  grown  by  plasma-assisted  MBE  in  a  modular  Gen  II  reactor  using  liquid 
gallium,  solid  Nd,  and  a  nitrogen  plasma.  The  photoluminescence  (PL)  spectrum 
from  the  lowest  4F3/2  multiplet  to  all  4L/2  multiplets  was  measured  by  a  free-space 
variable  amplifier  silicon  diode.  A  demonstration  of  the  coherence  properties 
requires  a  high-quality  sample  with  a  narrow  homogeneous  broadening.  For  this 
reason,  the  PL  linewidth  was  characterized  for  various  growth  conditions.  The 
experiments  were  carried  out  to  observe  coherent  effects,  and  various  experimental 
setups  are  described.  In  parallel  with  experiments  on  the  coherence  properties  of 
Nd3+:GaN,  we  also  investigated  the  effect  of  a  varying  electric  field  on  the  Nd  PL 
spectrum  to  provide  a  tunable  memory.  To  vary  the  applied  field,  we  designed  and 
grew  a  series  of  Nd-doped  GaN  p-i-n  structures,  strain-balanced  superlattice 
structures,  and  periodic  unbalanced  structures. 

2.  Introduction 


Because  of  an  increasing  demand  to  process  large  databases  in  real  time,  the  Army 
has  a  clear  interest  in  realizing  quantum  computing.  The  Army  also  has  a  great 
interest  in  new  approaches  to  secure  communications.  This  can  be  answered  by 
quantum  communication  schemes  using  polarization  states  of  photons.  Realization 
of  quantum  information  and  quantum  communication,  which  promise 
unprecedented  computational  capacities  via  quantum  computing  and  unconditional 
communication  security  via  quantum  cryptography,  stimulated  some  new 
approaches  to  optical  storage.1  Quantum  information  processing  requires  temporal 
storage  for  quantum  data.  The  coherent,  efficient,  and  reversible  mapping  between 
quantum  light  and  matter  represents  a  fundamental  challenge  in  the  field  of 
quantum  information  science.  Overcoming  this  challenge  would  enable  the 
realization  of  a  quantum  memory  for  light,  a  device  that  can  store  and,  on-demand, 
recall  quantum  states  of  light  with  high  efficiency  and  fidelity.  Such  memories  have 
potential  applications  for  use  in  a  quantum  repeater,  a  resource  allowing  for 
quantum  communication  over  long  distances.  Other  applications  include  linear 
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optics  quantum  computation,  deterministic  single  photon  sources  and  multiphoton 
quantum-state  engineering.  For  practical  applications,  it  is  likely  that  a  quantum 
memory  with  a  high  multimode  capacity  will  be  necessary. 

The  research  objective  of  this  Director’s  Research  Initiative  (DRI)  was  to  work  on 
developing  solid-state  quantum  memory  using  cryogenically  cooled  rare-earth-ion- 
doped  crystal,  Nd3+:GaN.  The  samples  were  grown  using  the  US  Army  Research 
Laboratory’s  (ARL’s)  molecular  beam  epitaxy  (MBE)  capabilities.  The  approach 
is  based  on  the  coherent  transient  interaction  between  light  and  the  ions  in  a 
semiconductor.  We  investigated  the  energy  level  structure  of  the  neodymium  (Nd) 
ions  embedded  into  the  semiconductor  crystal  lattice,  gallium  nitride  (GaN),  and 
developed  a  quantum  description  of  the  2-  and  3-pulse  photon  echo  effect.  To  do 
so,  we  identified  the  range  of  photon  frequencies  that  can  be  efficiently  stored  in 
the  quantum  memory.  Applied  electric  field  gradients  are  required  for  controlled 
inhomogeneous  broadening.  This  was  accomplished  using  semiconductor  p-n 
junctions  and  GaN/AlGaN  heterostructures,  which  allow  us  to  apply  much  stronger 
fields  than  currently  possible  in  quantum  memories  based  on  insulating  crystals  and 
glasses.  In  addition,  the  strong  internal  electric  fields  in  GaN  polar  heterostructures 
provide  a  unique  opportunity  to  control  the  frequency  of  the  electronic  transitions 
and  design  a  quantum  memory  with  a  multimode  capacity.  Also,  the  use  of 
semiconducting  host  materials  may  allow  optical  control  of  the  relevant  electronic 
transitions  and  field  control  through  screening  in  the  future  designs  of  quantum 
memories. 

3.  Approach  to  Achieving  Research  Objectives 

Many  different  approaches  toward  developing  quantum  memory  have  been 
proposed.  Recently,  significant  progress  has  been  achieved  using  atomic  ensembles 
for  manipulating  the  propagation  and  quantum  state  of  an  optical  field.2  The  rare- 
earth-ion-doped  crystals  (REICs)  are  promising  candidates  for  quantum  memory. 
REICs  are  characterized  by  large  optical  inhomogeneous  broadening  that  enables 
storage  and  recall  of  coherent  information  by  manipulating  and  controlling  the 
inhomogeneous  dephasing  using  echo  techniques.  Although  traditional  photon 
echoes  are  not  directly  used  in  the  single-photon  regime  because  of  spontaneous 
emission  noise  induced  by  the  Ti-pulse,3  photon  echo  techniques  avoiding  this  noise 
have  been  proposed:  controlled  reversible  inhomogeneous  broadening  (CRIB)4  and 
more  recently  atomic  frequency  combs  (AFCs).5 
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3.1  Quantum  Memory  Based  on  2-Pulse  Photon  Echo 


As  we  mentioned  previously,  the  traditional  photon-echo-based  memory  generally 
suffers  from  strong  limitations  when  used  for  storage  of  single  photons;  therefore, 
it  cannot  be  directly  used  for  quantum  communication.6  However,  the  standard 
photon  echo  studies  are  a  critical  first  step  and  give  rise  to  the  development  of  the 
quantum  storage  approach  based  on  CRIB,  which  is  well  suited  for  storage  of  single 
photons.  We  therefore  give  a  brief  description  of  the  traditional  2-pulse  photon  echo 
technique. 

The  underlying  principle  of  quantum  memory  based  on  the  photon  echo  technique 
can  readily  be  explained  in  a  Bloch  vector  representation  since  it  allows  for  a 
geometrical  interpretation  of  the  density  matrix  evolution.  Using  Pauli  matrix 
decomposition,  one  can  map  the  evolution  of  the  2-level  system  (representing 
selected  states  of  the  rare-earth  ion  in  the  solids)  into  rotation  of  the  3-D  Bloch 
vector  B  =  (u,v,w)  =  (pl2+p2l,Kp2l-pl2),p22-pn),  where  p  11,22  is  the  population  of 

the  states,  and  pn  is  the  coherence.  Figure  1  shows  a  schematic  of  the  2-pulse 
photon  echo  in  the  Bloch  vector  representation. 


Fig.  1  Evolution  of  the  Bloch  vector  during  the  2-pulse  photon-echo  process,  a)  n/2  pulse 
rotates  the  Bloch  vector  from  the  negative  w  axis  along  the  u  direction,  b)  The  individual  Bloch 
vectors  precess  in  the  uv  plane  and  dephase  due  to  the  inhomogeneous  broadening  of  the 
transition,  c)  At  time  t  =  t,  a  71-pulse  rotates  all  Bloch  vectors  around  the  v  axis,  d)  The  Bloch 
vectors  realign  and  build  up  a  macroscopic  coherence  that  results  in  a  photon-echo  emission 
at  time  t  =  2r. 

In  some  cases,  analytic  calculations  of  the  echo  signal  intensity  can  be  carried  out 
using  coupled  Maxwell-Bloch  equations.7  Taking  into  account  the  inhomogeneous 
distribution  of  the  transition  frequency,  one  can  find  the  intensity  of  the  echo  signal 
P(t)  oc  sin  (A1)sin4(A!  /2)e"2r'e_<72(,“2r)2/2,  where  A*  is  the  pulse  area  of  the  Upulse,  k 

=  1,2,  y  is  the  dephasing  rate,  and  a  is  the  width  of  the  inhomogeneous  distribution. 
The  dimensionless  pulse  area  is  defined  as  a  time  integral  over  the  electric  field 
amplitude  Eo(t)  multiplied  by  the  electric  dipole  moment  of  the  interlevel  transition 
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di2  and  divided  by  Plank’s  constant.  It  is  not  necessary  that  the  pulse  areas  be  equal 
to  7t/2  and  n,  although  these  areas  lead  to  a  maximal  signal.  What  is  necessary  is 
that  the  second  pulse  produces  at  least  a  partial  reflection  about  the  uw  plane.  This 
reflection  takes  the  Bloch  vector  components  u  +  iv  into  u  -  iv  or,  equivalently, 
takes  density  matrix  element  pn  into  pu.  Since  pn  and  pi\  are  related  to  the  real 
and  imaginary  parts  of  the  average  dipole  moment  operator,  the  second  pulse  must 
couple  these  real  and  imaginary  parts.  Such  coupling  is  impossible  for  a  linear 
atom-field  interaction.  Thus,  by  its  very  nature,  the  photon  echo  can  occur  only 
when  a  nonlinear  atomic  response  is  present. 

3.2  Quantum  Memory  Based  on  Controlled  Reversible 
Inhomogeneous  Broadening  (CRIB) 

A  few  decades  ago,  researchers  demonstrated  how  to  benefit  from  the  large  ratio 
between  the  homogeneous  and  inhomogeneous  linewidths.8  The  CRIB  photon- 
echo  technique  is  a  4-wave  mixing  technique,  with  pulses  being  separated  in  time. 
The  3  input  pulses  manipulate  the  system,  putting  the  absorbers  in  the  material  in 
superposition  states,  which  at  some  given  time  all  oscillate  in  phase  creating  a 
macroscopic  dipole  moment.  The  output  pulse  (echo)  is  emitted  by  this  temporally 
created  macroscopic  dipole  moment.  To  store  a  sequence  of  data,  a  short 
preparation  (write)  pulse,  Ew(t),  with  a  pulse  area  equal  to  tc/2  and  a  bandwidth 
equal  or  larger  than  the  required  data  rate  is  sent  into  the  sample  (Fig.  2).  The  data 
pulse  sequence,  Edit),  at  a  maximum  data  rate  equal  to  the  bandwidth  of  the  pulse 
Ew(t),  and  duration  shorter  than  the  homogeneous  relaxation  time  of  the  absorbers, 
is  then  sent  into  the  media.  The  input  radiation  (both  write  and  data  pulses)  interacts 
coherently  with  the  absorbers.  As  a  result,  the  frequency  spectrum  of  the  input  data 
sequence  (amplitude  as  well  as  phase)  becomes  imprinted  as  a  frequency-dependent 
modulation  onto  the  absorption  profile.  Since  both  amplitude  and  phase  are  stored, 
the  full  temporal  information  of  the  light  pulses  is  contained  in  the  spectral 
interference  patterns  engraved  in  the  sample,  and  later  pulses  can  interact  with  these 
spectral  gratings  to  perform  read-out.  Mathematically,  it  can  be  expressed  by  the 
following  Fourier  transform  equation.  When  the  electromagnetic  fields  of  the  3 
input  pulse  sequences,  Ew(t),  Ed(t),  and  Er(t),  are  far  from  saturating  the  optical 
transition  and  the  time  separation  between  the  pulse  sequences  is  smaller  than  the 
transition  relaxation  time,  the  output  field,  Ee(t),  as  a  function  of  time  is 
Ee  ( t )  oc  J  E*w[a>)  Ed  (&>)  Er  ( co)emda ).  It  is  clear  that  the  echo  field  is  just  a  copy  of  the 

Ed(t)  when  the  spectra  of  Ew(co)  and  Er(co)  are  constant  across  the  data  bandwidth, 
Ed(  co). 
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Fig.  2  a)  Stimulated  photon-echo  schematic,  b)  Generalization  of  the  3-pulse  echo  to  data 
storage. 


The  CRIB  approach  involves  broadening  an  initially  narrow  absorption  line  using 
the  linear  Stark  effect  and  an  applied  external  electric  field  gradient.9  The  width  of 
the  broadened  line  should  match  the  spectral  width  of  the  light  that  is  to  be 
absorbed.  Read  out  or  re-emission  of  the  stored  light  is  triggered  by  changing  the 
sign  of  the  external  electric  field,  which  inverts  the  atomic  transition  frequencies 
around  the  central  frequency.  The  rephasing  mechanism  can  be  understood  by  the 
following  argument.  Atoms  at  detuning  A  with  respect  to  the  carrier  frequency  of 
the  light  accumulate  a  phase  -zAx  between  absorption  at  t  =  0  and  the  polarity 
inversion  of  the  electric  field  at  t  =  x.  By  inverting  the  polarity,  the  atoms  invert 
their  detunings  A  — ►  -A.  After  a  time  t  =  2r  the  atoms  will  have  accumulated  an 
opposite  phase  factor  zAx  that  results  in  a  rephasing  of  all  atomic  dipoles.1  A 
schematic  of  the  CRIB  protocol  is  shown  in  Fig.  3. 


Step  1  Prepare  narrow  absorption  line 


Step  2:  Broaden  line 


A 


3 

11  1  1 

11  (l 

p  ' 

In 
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< 

HH 
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— =► 

Linear  Stark  shifts  with  external  electric  field  ® 
Step  3:  Absorb  a  photon 

Step  4:  Trigger  re-«ni$sion 


A 


CO 

CRIB  by  diaugiiis!  the  polarity  of  the  external  electric  field 


Fig.  3  Quantum  state  storage  based  on  CRIB  in  crystals  with  an  ensemble  of  ions  having  a 
broad,  inhomogeneously  broadened  absorption  line 
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4.  Electronic  States  in  Rare  Earth  Ions 


In  this  study  we  use  Nd3+:GaN.10  Rare-earth-ion-doped  solids  are  of  particular 
interest  because  of  the  absence  of  motion  of  the  absorbers,  as  well  as  their  excellent 
optical  and  hyperfine  coherence  times  at  cryogenic  temperatures,  which  promise 
long  storage  times.  As  a  result  of  different  strains  in  the  host  medium,  they  also 
exhibit  large  inhomogeneous  broadening  on  their  optical  transitions,  which 
provides  the  fast  decay  of  excited  atomic  coherence  upon  absorption  of  broadband 
light.  The  ratio  between  the  inhomogeneous  and  homogeneous  transition  linewidths 
can  get  up  to  108  in  REICs.11  Therefore,  by  using  the  frequency  dimension  to 
address  and  store  optical  data,  one  could  potentially  store  and  address  thousands  of 
data  bits  at  a  single  spatial  location. 

Material  requirements  for  CRIB-based  quantum  memory  include  the  possibility  to 
inhomogeneously  broaden  an  optical  transition  in  a  controlled  way.  The  induced 
broadening  needs  to  be  large  compared  to  the  inherent  broadening  of  the  transition, 
and  it  should  be  possible  to  reverse  it  in  a  short  amount  of  time  compared  to  the 
inverse  of  the  inherent  broadening.  The  shielding  of  the  4/  electrons  reduces 
electron  orbit-lattice  interaction  as  compared  to  other  solids,  thereby  minimizing 
the  effect  of  dynamic  perturbations  by  phonons.  This  results  in  a  transition  intensity 
concentrated  in  narrow  zero-phonon  lines.  For  crystalline  hosts  the  homogeneous 
linewidths  are  typically  around  a  few  kilohertz.1 

The  ion  energy  levels  can  be  manipulated  using  the  Zeeman  and  Stark  effects.  Of 
particular  interest  for  CRIB  is  the  linear  Stark  effect.  For  an  impurity  ion  in  an 
inversionless  crystal  field,  the  states  of  different  parity  are  mixed  and  thus  a  linear 
Stark  effect  in  an  external  electric  field  is  possible.12  Provided  the  permanent 
electric  dipole  moments  of  the  states  coupled  by  the  optical  field  are  different,  it 
leads  to  a  shift  in  resonance  frequency,  which  can  be  exploited  for  controlled 
reversible  inhomogeneous  broadening  of  the  associated  transition.  The  magnitude 
of  the  Stark  shift  depends  on  the  particular  transition,  the  crystal  host,  and  the 
orientation  of  the  electric  field  with  respect  to  the  permanent  dipole  moment 
difference.  The  largest  Stark  shifts  are  of  the  order  of  100  kHz/Vcm-1.12  In  RE- 
doped  crystals,  where  the  application  of  a  DC  electric  field  leads  to  a  shift  or 
discrete  splitting  of  an  absorption  line,  controlled  inhomogeneous  broadening  on 
the  order  of  several  hundred  megahertz  can  potentially  be  induced  by  applying  an 
electric  field  gradient  with  a  field  strength  varying  between  -1000  VcirT1  and 
+1000  Vcnr1. 
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Nd3+  has  three  4/ electrons  (14  electrons  are  needed  to  completely  fill  the  shell). 
Labeling  the  various  levels  is  based  on  the  Russell-Saunders  (RS)  coupling  scheme. 
In  the  RS  scheme  the  interaction  between  orbital  and  spin  angular  momenta  of  the 
electrons  is  neglected  in  the  first  approximation,  and  the  total  orbital  angular 
momentum  is  quantized. 

The  orbital-momentum  vectors  of  2  electrons  l\  and  h  give  the  total  momentum  L 
with  values  from  I1+L2  to  I/1-/2I.  The  4/ electron  has  1  =  3.  Each  value  of  L  is  similarly 
coupled  to  h,  and  so  on.  The  total  orbital  momenta  L  =  0,1,2,3,4,5,6,7,8,9,10,11 
corresponds  to  electronic  terms  S,P,D,F,G,H,I,K,L,M,N,0,Q.  Similarly,  the  spin 
momenta  of  the  individual  4/ electrons  are  vectorially  added  to  give  S. 

5.  Review  of  Previous  Growth  of  Neodymium  (Nd)-Doped 
Gallium  Nitride  (GaN) 

Nd  has  previously  been  incorporated  into  GaN  through  ion  implantation  and 
reactive  cosputtering  deposition.  Ion  implantation  has  limitations  due  to  damage 
effects  and  doping  profiles.  In  many  cases,  a  postimplantation  anneal  is  required 
because  of  the  amorphization  of  the  GaN  lattice;  however,  annealing  does  not 
completely  recover  the  damage.  For  ion-implanted  Nd  in  GaN,  the  majority  of  Nd3+ 
ions  sit  on  the  Ga  substitutional  “site”;  however,  4  other  sites  associated  with 
impurities  and  defects  are  also  identified  within  the  GaN  matrix.  For  reactive 
cosputtering  deposition  of  GaN  with  RE  atoms,  a  postgrowth  anneal  (>700  °C)  in 
nitrogen  is  also  necessary  for  photoluminescence  (PL)  emission.  As  part  of 
previous  DRI  work  at  ARL,  we  successfully  achieved  in  situ  Nd  doping  of  high 
crystal  quality  GaN  grown  by  plasma-assisted  MBE.13  The  GaN:Nd  samples, 
without  any  postgrowth  annealing,  exhibited  strong  PL  emission  under  direct 
pumping  of  the  Nd  centers  that  could  be  identified  with  the  4F3/2  ->-4F/2,  4Fm 
— »4/n/2,  and  4Fyi  — >4/i 3/2  multiplets  (Fig.  4).  Moreover,  the  scarcity  and  relative 
weakness  of  the  emission  from  extra  PL  peaks  from  above  bandgap  pumping 
suggest  improved  incorporation  at  the  Ga  substitutional  sites  and  allowed  us  to 
resolve  the  Stark  energy  sublevels  of  the  4/ states  using  PL  and  photoluminescence 
excitation  (PLE)  spectroscopies.10 
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Fig.  4  Low-temperature  (~30  K)  PL  from  900  to  1500  nm  for  GaN:Nd  pumped  at  836  nm 
showing  the  characteristic  Nd  transitions.  (Reprinted  from  Readinger  et  al.13) 

The  GaN  layers  were  grown  on  single-side  polished  opiane  sapphire  by  MBE  using 
standard  effusion  cells  and  a  UNI-Bulb  RF  plasma  source  for  activated  nitrogen. 
Following  the  sapphire  nitridation,  an  aluminum  nitride  (AIN)  buffer  layer  of 
approximately  25  nm  was  grown.  The  buffer  was  followed  by  a  200-nm  undoped 
GaN  base  layer  and  a  Nd-doped  approximately  1-pm  GaN  layer  with  various  Nd 
concentrations  up  to  approximately  5  atomic-percent  (at%).  The  Nd  cell 
temperature  was  adjusted  while  maintaining  fixed  growth  conditions  (growth 
temperature,  Ga  flux,  nitrogen  flux,  rf  power,  nitrogen  gas  flow  rate,  etc.)  to  control 
the  dopant  incorporation. 

The  material  grown  for  this  work  was  analyzed  using  Rutherford  backscattering 
spectroscopy  and  secondary  ion  mass  spectrometry  (SIMS)  for  quantitative 
analysis  of  the  Nd  concentration,  X-ray  diffraction  (XRD)  for  structural  quality, 
and  PL,  PLE,  reflectance,  and  transmittance. 

6.  Theoretical  Analyses  of  2-  and  3-Level  Models  and  Design 
Parameters 


6.1  2-Level  System 


In  the  rotating  wave  approximation  (RWA),  the  Hamiltonian  has  the  form 


H  =  - 


h 

2 


r0 

VQ 


Q  " 
2A<yy 


(1) 
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where  A  G)  -  (0-  E2  /  tl  is  the  detuning  between  the  transition  frequency  and  the 
center  frequency  of  the  laser  field,  taken  as 
E(t )  =  E0(t)cos(o)t  +  ^0),  Q  =  E0(t)dn  /  tl  is  the  Rabi  frequency,  du  is  the  dipole 
moment.  The  model  is  illustrated  in  Fig.  5. 


Aoj  !x 


1 2} 


fl(0 


|l) 


Fig.  5  Schematic  of  a  2-level  system 

To  take  into  account  spontaneous  decay  of  the  exited  state,  it  is  constructive  to 
rewrite  the  density  matrix  equation  (/HO  =  ~  P~  pH] )  in  the  Liouvillian 
form 


X(t)  =  -(4  +  L)X(t), 


(2) 


where 


"a2(0" 

X(t)  = 

Pi  i  (0 
Ai(0 

\  Pl2  (0  y 

the  dynamical  part  of  the  Liouvillian  is 

( 
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(3) 


(4) 


the  relaxation  and  dephasing  are  taken  into  account  by  using  the  superoperator 
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0  0  0 


0—0  0 


0  0  0  - 

T 


0  0  0  — 
T, 


i/r2*  =i/(27;)+r  is  the  dephasing  rate,  and  T\  is  the  life  time  of  the  excited  state. 

The  exact  form  of  the  dynamic  equation  for  the  density  matrix  elements  is 

P\ 2  P\ 2  (  1/^2  )  (/^ll  Pll  )  ^  ^ 

P21  —  P21  (i^G) -1/T2 )  +  (/?n  —  P22 ) ^ ^ 2, 

(6) 

Al  —  P'2'2  +  (/^21  _  P\2  )  ^ 

P22  ~  ~~ P22  ^  T\  ~  (/?21  ~  P\2  )  ^  ^ 

Taking  into  account  that  /?n(0  +  /^(O  =  1’  we  find  the  steady-state  solution 


Pi  2 
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n 
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Q 

A®+i/r2* 
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A  co2 
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Q  %/T* 

2A<a2+i/r2‘2+Q2r1/r2* 
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cPtjt* 

~2 

A<y2 

+  1/T*2  +n2Tl/T*' 

/  ’  .  — 


We  see  that  in  the  resonant  case,  =  0,  the  population  of  the  excited  state  is 

1  Q  2TT* 

PS2 2= - (8) 

2i+Q27;r2 

and  in  the  limit  Q2TiT2*  »  1  the  state  populations  of  the  2-level  system 

Pu  ~  P22  =1/2. 
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6.2  Bloch  Equations 


Using  pn  =  p2\  -  ( u-iv )  /  2,  p22  —pn  =  n,  and  p12  +  pu  - 1,  we  can  rewrite  the 
density  matrix  (Eq.  6)  in  the  Bloch  equation  form  as 

u  =  -u/T2  +  Aryv, 

v  =  -v  /  T2  -  Acou-Qn,  (9) 

n  =  (n  + 1)  /  Tj  +  Qv. 


6.3  3-Level  System  and  Fast  Nonradiative  Decay  of  the  Top 
Level 


To  describe  fluorescence  from  an  additional  level  I3>  populated  by  fast 
nonradiative  transition  from  the  excited  state  I2>,  we  consider  the  3-level  system 
shown  in  Fig.  6. 


Fig.  6  Schematic  of  a  3-level  system  with  nonradiative  decay  of  the  top  level 


In  this  case,  we  have  to  include  an  additional  equation  of  motion  of  the  population 
of  the  excited  state  I3>.  The  dynamical  part  of  the  Liouvillian,  Eq.  4,  does  not 
change  while  the  relaxation  part,  Eq.  5,  becomes 


L  = 


UT1+ynl2 
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(10) 


where  /„  is  the  nonradiative  decay  rate  of  the  state  I2>,  and  7Y3'  is  the  lifetime  of 
the  state  I3>. 

The  system  of  the  dynamic  equation  for  the  density  matrix  elements  becomes 
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Pn  =  Pn  (~iAco ~ 1 1  Ti  ~  A  /  2)  -  (  Ai  -  P22 ) in  1 2 , 

Ai  =  Ai ' 1  /  A*  -  A  /  2)  +  (a  1  -  A2 )  /  2, 

Ai  =  Pn 1 T1  +  As  / 7’i<3)  +  ( Ai  -  A2 ) /  2,  (1 1) 

A:  =  —  P22  (1  /  +  A  )  —  (  Ai  —  A2 )  (  2, 

A3  =  A2A  —  A3  '• 


»  JJ  I  LL 


Taking  into  account 

,  = _ /7;iQ(^+l)(2  +  7;:t(ri,+2/A^)) _ 

A2  n2TX(rX3) + 2K/X  +  2) + (a*; + D(4A«2r;2 + for; + 2)2) 

.  = _ +  l)(2t  +  A*  O' A  +  2Ato)) _ 

A2  n^r;  (r„r,(3)  +  2)(r„r; + 2) + for,  +  i)(4A«2r;2  +  for; + 2)2 )  ’ 

,  (a?! +i)(4A^2r;2 +(7„r; +2)2)+Q2r1r;(7j; +2) 

Q^for® + 2)(rX  +  2) + for,  +  i)(4A®2r2*2 + for;  +  2)2)  ’ 

= _ tfrtiirX+Q _ 

A2  QV;  for/3)  +  2)for;  +  2) + for,  +  i)(4Afor2*2  +  for;  +  2)2 )  ’ 

yjXTi^Wi  +  2) 

n2r,r2* for,(3) + 2)(a A* + 2) + fo^i + i)(4Afor;2  +  for;  +  2)2) ' 
Note  that  Eq.  12  is  identical  to  Eq.  7  in  the  limit  Jn  A  0. 


(12) 


In  the  steady-state  regime,  from  the  equation  for  p33,  Eq.  11,  we  see  that 

P33  =  * 3 1 Y ,1 P22  •  Therefore,  now  we  consider  steady-state  expression  for  the  p'l2- 

For  the  resonant  excitation,  Aco  =  0,  we  obtain 


Q2772* 


Pl 2  ( rX3)  +  2) + (a7!  + 1)(  A  A* + 2) ' 

Assuming  that  ynTi  »  1  and  7]<3)  «  7j  ,  we  find 
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(13) 
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(15) 
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6.4  3-Level  System  and  Fast  Nonradiative  Decay  of  the 
Intermediate  Level 

Figure  7  shows  a  schematic  of  the  3-level  system  with  fast  nonradiative  decay  of 
the  intermediate  level  I3>. 


Fig.  7  Schematic  of  a  3-level  system  with  nonradiative  decay  of  the  intermediate  level 

In  this  case,  the  dynamical  part  of  the  Liouvillian,  Eq.  4,  does  not  change  while  the 
relaxation  part,  Eq.  10,  becomes 
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where  1  /  T23  is  the  radiative  decay  rate  of  the  state  I2>  to  the  state  I3>,  and  Yn  is 
the  nonradiative  decay  rate  of  the  state  I3>.  Now  1  /  T2  =  1  /  2(1  /  T{  +  l/T23)  +  T. 

The  system  of  the  dynamic  equation  for  the  density  matrix  elements  becomes 

Pi  2  —  Pi  2  O'AfiJ  —  1/7;  )  —  (/?!!  —  P22  )  /  2, 

p2\  —  P2\ i^co -HT2  )  +  (/?jj  —  P22 ) / 2, 

Pu-Pi2ITx  +  PiP  11  +  ( Pi  1  —  Pi  2 )  ^2  /  2,  (17) 

P22  —  ~ P22  A  ^  ^1  +  1  /  7)3  )  —  (/+!  —  /?!  2  )  /  2, 

7*33  =  P22  !  ^23  —  P33Y n  * 

Again,  taking  into  account  that  /?33  +  A2  Ai  =  ^  we  find  the  steady-state  solution 
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Pn  = 


/r;Q(r„r1+i)(2+r;(r„  +  2/A®)) 
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(18) 


In  the  steady-state  regime,  from  the  equation  for  <933,  Eq.  17,  we  see  that  gs 22  = 
T23yngs33.  Therefore,  now  we  consider  the  steady-state  expression  for  osn. 

For  the  resonant  excitation,  A co  =  0,  we  obtain 
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Q->'2 


7„r,j?r23n2 


27n(r,  +  r23)  +  r,7?n2(27„r23  + 1)  ■ 


(19) 


Assuming  that  ynTi3  »  1,  we  obtain 


1  T^O2 

2T1/r23  +  i  +  r1r^2- 


(20) 


The  steady-state  population  of  the  excited  level  is  very  similar  to  the  2-level  system 
case,  Eq.  8,  in  the  limit  Q2T172*— »  00  the  state  populations  gs\  \  =  gs 22  =  1/2.  The 
difference  is  that  now  we  have  population  inversion  on  the  I2>- 13>  transition,  since 
gs 33— »  0;  to  be  exact,  it  goes  as  1  /T i3yn,  which  is  much  less  than  1. 

7.  Steady-State  Population 

Figure  8  shows  the  steady-state  population  in  Eq.  8  as  a  function  of  the  Rabi 
frequency.  As  mentioned  previously,  the  population  of  the  states  is  equal  to  1/2  at 
the  limit  Q  XiT*2  —XX). 


Approved  for  public  release;  distribution  unlimited. 


14 


Introducing  a  new  variable  x  =  QxTiT*2,  we  see  that  the  steady-state  population  in 
Eq.  8  is  f{x)  =  0.5x2/(l+x2).  The  second  derivative  (Fig.  9)  of  this  function  is  equal 
to  zero  at  xo  =  1/V3  (therefore,  QO  =  \N?>T\T*i),  which  can  be  used  to  determine 
the  dipole  moment  from  the  experimental  data.  Figure  10  shows  the  steady-state 
population  in  Eq.  15  as  function  of  the  Rabi  frequency.  We  assume  that  T\  ~  Tiy)\, 
YT\  »  1,  and  yJT  =  10.  For  comparison,  the  dashed  line  shows  the  steady-state 
population  in  Eq.  8. 


n^fr i  n 


Fig.  8  Steady-state  population  in  2-level  system,  Eq.  8 


Fig.  9  Useful  dependencies  to  analyze  the  steady-state  population  as  a  function  of  the  field 
amplitude 
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Fig.  10  Steady-state  population,  see  Eqs.  8  and  15 

Figure  1 1  shows  fix)  =  x2/{a2  +  x2)  (a2  =  ynT *2  +  2),  its  first  and  second  derivatives 
as  a  function  of  x  =  qVt,3)i  T* 2  .  The  second  derivative  of  this  function  is  equal  to 
zero  at  xo  =  a/V3;  therefore,  Qo  =^lynT* 2  +  2/^3T\ T*i. 


X 


Fig.  11  Useful  dependencies  to  analyze  the  steady-state  population  as  a  function  of  the  Held 
amplitude 

Figure  12  shows  the  steady- state  population  in  Eq.  20  as  a  function  of  the  Rabi 
frequency. 

We  assume  that  T\  ~  T23,  and  ynT22 »  1.  For  comparison,  the  dashed  line  shows  the 
steady-state  population  in  Eq.  8.  Figure  13  shows  /(x)  =  0.5x2/(a2+x2)  ( a 2  = 
T\/T 23+I),  its  first  and  second  derivatives  as  a  function  of  x  =  Q^T\T*2 .  The  second 
derivative  of  this  function  is  equal  to  zero  at  xo  =  a/V 3;  therefore,  Qo  =  ^T\/T2i  + 

1/V3  7)7^2. 
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n^fr i  rf 


Fig.  12  Steady-state  population,  see  Eqs.  8  and  20 


Fig.  13  Useful  dependencies  to  analyze  the  steady-state  population  as  a  function  of  the  field 
amplitude 


8.  Pulsed  Excitation  of  the  2-Level  System 


For  the  2-level  system  wave  function,  lvF(t)/)  =  a\W  /+ci2l2),  in  the  Schrodinger 
picture  the  Hamiltonian  has  the  form 


h  (  0  + 


(21) 


where  (02  =  E?/h  is  the  transition  frequency  (we  chose  the  energy  of  the  ground  state 
Ei  =  0),  co  is  the  center  frequency  of  the  laser  field,  Q.(t)  =  E(t)d\2/h  is  the  time- 
dependent  Rabi  frequency,  which  is  determined  by  the  laser  field  pulse  envelope, 
E(t),  and  the  dipole  moment,  du.  Using  interaction  representation  and  RWA,  we 
obtain  the  Hamiltonian  as 
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H  =  -p(t) 


0 


(22) 


where  detuning  Aco  =  co-coi. 

To  take  into  account  relaxation,  we  use  the  master  equation  for  the  density  matrix 

pm-- 

s(t  )  =  -ji\He-eH\).  (23) 

which  allows  us  to  correctly  include  relaxation  terms. 

In  the  Liouvillian  form,  the  final  master  equation  has  the  form 

X(t)  =  L{t)X(t)  .  (24) 


where 


x(t)  = 


^  Oi2(0  ^ 

021(f) 

-11  (f ) 

92 2(f)  j 


(25) 


—iAuj  —  J- 


m  = 


0 


2  — T~ 


0 

mi 

2 

m> 


iAu-A.  m±  _Mu 

44-lkiL-  lTV  ij  13 


2 

tfl(t) 


Tj 


0  -J- 


"l\  / 


(26) 


First,  let  us  consider  the  limiting  case  of  no  relaxation  and  constant  Rabi  frequency 
Q.  With  the  initial  condition 
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(27) 


(  eu(0)  =  o  \ 


X(0)  = 


021 (0)  =  0 
011(0)  - 1 


022(0)  =  0  / 


we  obtain 


X(t)  = 


f  Q{ Aw(ccs  tA- 1 ) -WAul+n?  sin  tA  )  ^ 
2fA^-  -£J-L 

H  (  All.'(cos  t A  —  1 )  -Hv* Aio^-t-fi2  si  n  t A  J 

2(A^-  -  £>-:. 


2Aw3+Q-+fi2  cos  tA 
2(A(Jj2+fia) 

nW3(£) 

Aufi-W  ' 


/ 


(28) 


Therefore,  in  the  resonant  case,  A co  =  0,  the  population  of  the  excited  state  is  Quit) 
=  sin2(tfi/2).  Since  the  PL  is  proportional  to  the  excited-state  population,  the 
integrated  PL  is 


IpL(t)  oc  jf  sin3  (f'ft/2)  dtf  =  ^  j  .  (29) 

Figure  14  shows  the  excited-state  population,  Eq.  28,  and  the  integrated  PL,  Eq.  29, 
as  a  function  of  time. 


t 


Fig.  14  Integrated  PL  and  the  excited-state  population  vs.  time  for  the  limiting  case  of 
constant  Rabi  frequency 
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To  show  the  relaxation  effect  and  take  into  account  the  pulse  envelope,  we  solved 
the  master  equation  numerically.  Figure  15  shows  the  excited-state  population  (left 
panel)  and  the  integrated  PL  (right  panel)  as  a  function  of  time  in  the  case  of  pulse 
excitation  with  the  Rabi  frequency  envelope  shown  in  Fig.  16. 


0  2  4  6  8  10  12  14  0  2  4  6  8  10  12  14 

t  t 

Fig.  15  Integrated  PL  and  the  excited-state  population  vs.  time  for  the  pulse  excitation 


Fig.  16  Rabi  frequency  vs.  time  for  the  pulse  excitation 


9.  Experimental  Work  for  Nd3+:GaN  Director's  Research 
Initiative  (DRI) 


9.1  Growth  of  Nd3+:GaN 

Epitaxial  GaN:Nd  was  grown  by  plasma-assisted  molecular  beam  epitaxy  (PA- 
MBE)  in  a  modular  Gen  II  reactor  using  liquid  gallium,  solid  Nd,  and  a  nitrogen 
plasma.  A  gallium  beam  equivalent  pressure  of  1  x  10-7  Torr  and  a  N2  gas  flow  of 
0.40  seem  was  used  to  grow  stoichiometric  GaN  at  180  nm/h  with  a  step-flow 
growth  morphology  resulting  in  a  surface  roughness  of  less  than  1  nm.  A  series  of 
Nd  doping  calibrations  were  grown  that  consisted  of  3  layers  of  GaN:Nd  with 
different  Nd  source  temperatures.  The  resulting  Nd  concentrations  ranged  between 
5  x  1017  to  5  x  1018  cm-3.  Once  calibrated,  GaN:Nd  films  were  grown  on  free¬ 
standing  GaN  substrates  with  a  target  Nd  concentration  of  approximately  0.01%, 
or  4  x  1018  cm-3. 
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9.2  Continuous-Wave  (CW)  Steady-State  Saturation 
Photoluminescence  (PL)  Setup  and  Experiment 


Previous  experimental  work  to  characterize  the  emission  and  excitation  spectrum 
of  Nd3+  ions  doped  in  GaN  has  been  used  to  identify  transitions  between  the  4F3/2 
excited  state  to  the  4Ig/2,  4Ii  1/2,  and  4Ii3/2  multiplets  (Fig.  17).  For  the  purpose  of 
quantum  memory  applications,  it  was  necessary  to  identify  a  2-level  system  for 
study  within  the  rich  Nd  spectrum. 


Fig.  17  Simplified  energy  diagram  of  the  manifolds  in  the  Nd3+  ion  in  GaN.  (Reproduced 
from:  Metcalfe  G,  Redinger  E,  Shen  H,  Woodward  N,  Dieroff  V,  Wraback  M.  Crystal-field 
split  levels  of  Nd3+  ions  in  GaN  measured  by  luminescence  spectroscopy.  J  Appl  Phys. 
2009;  105:053101.,  with  the  permission  of  AIP  Publishing.) 
doi:  http://dx.doi.Org/10.1063/l.3082500 

To  characterize  Nd3+:GaN  samples,  a  near-infrared  PL  setup  was  constructed  to 
study  the  luminescence  spectra  in  an  ARL-grown  1.2- pm- thick  Nd3+:GaN  sample 
with  0.51%  Nd  concentration  grown  by  PA-MBE  using  CW  and  time-resolved 
techniques  (Fig.  18,  left  panel).  Samples  were  excited  using  a  5-W,  532-nm 
continuous  wave  (CW)  laser  to  pump  a  CW  Ti:Sapphire  laser  tunable  between  700 
and  1000  nm.  A  CW  laser  system  was  used,  as  opposed  to  a  tunable  ultrafast  laser 
system  possessing  broadband  emission,  because  of  the  available  narrow  linewidth 
(<40  GHz),  which  enables  the  individual  transitions  to  be  excited  selectively.  The 
laser  is  passed  through  an  optional  electro-optic  modulator  (EOM)  or  acousto-optic 
modulator  (AOM),  which  can  modulate  a  square-wave  temporal  profile  on  the 
incident  CW  laser  beam  up  to  20  MHz.  The  laser  is  incident  on  the  sample,  which 
is  mounted  to  a  cold  finger  inside  of  a  closed-cycle  helium  (He)  cryostat  capable  of 
providing  sample  temperatures  between  10  K  and  300  K.  All  measurements  for  this 
study  have  been  made  at  the  10-K  low-temperature  limit.  The  PL  emission  is  then 
collected  using  a  pair  of  lenses  to  collimate  and  focus  the  collected  light  onto  the 
entrance  slit  of  an  850-mm  spectrometer  with  a  600-g/mm  grating  and  600-pm- 
wide  entrance  and  exit  slits.  After  the  exit  slit  of  the  spectrometer,  a  rotatable 


Approved  for  public  release;  distribution  unlimited. 

21 


concave  mirror  allows  for  several  detection  options  of  the  PL  signal.  To  change 
between  multiple  detectors  with  a  range  of  sensitivity  and  bandwidth,  we  have  fiber 
coupled  one  of  the  outputs  so  that  we  can  easily  use  several  fiber-based  detector 
options  interchangeably.  A  multimode  fiber  allows  for  70%  coupling  efficiency 
from  the  spectrometer.  For  the  CW  PL  measurements,  both  a  sensitive,  amplified 
silicon  detector  and  a  liquid-nitrogen-cooled  Ge  detector  were  used  to  obtain  the 
PL  spectrum  using  a  lock-in  amplifier  synchronized  to  the  modulation  of  the 
excitation  laser. 


■-i 

i 

j 

Fiber  based  high  speed  51 
or  InGafo  PMT  detector 

J. 

Fiber  II 

1 

i 

wricun  1 

■ - 

Monochromator 

GOQ’g/mm 

j  i . fl  " 

He  Crvostst  ■  !  i - 

-It*  ■  g-'A  ;  |  COM  H — 


5W 
532  nn 


TiLSqpphire 
700-LOOO  nm  400  m W 


Fig.  18  Left  panel:  Experimental  setup  for  CW  PL  and  time-gated  PL.  Right  panel: 
Emission  spectrum  from  the  lower  4F3/2  state  to  the  4I<>/2  ground  multiplet;  the  arrow  indicates 
the  transition  between  the  ground  levels  of  each  state  at  916.2  nm,  which  is  the  strongest. 

The  PL  spectrum  from  the  lowest  4F3/2  multiplet  to  all  %/2  multiplets  was  measured 
by  a  free-space  variable  amplifier  Si  diode  and  is  shown  in  the  right  panel  of  Fig.  18. 
To  avoid  interference  from  laser  scattering,  we  initially  pumped  the  Nd3+:GaN 
sample  into  the  4Fs/2  state  and  looked  at  the  emission  from  the  4F3/2  state.  The 
integrated  emission  from  the  4F3/2  excited  state  to  the  %/2  ground  state  is  plotted  as 
a  function  of  the  exciting  electric  field  (calculated  from  the  incident  laser  intensity) 
in  the  left  panel  of  Fig.  19,  which  has  been  shown  capable  of  determining  the  dipole 
moment  of  the  transition  from  theoretical  calculations.  The  PL  emission  intensity 
is  fitted  assuming  a  functional  form  of  the  saturation  dependence  of 
Y=A(Bx2)/(1+Bx2).  The  resulting  fit  is  plotted  as  the  solid  blue  line  in  Fig.  19  (left 
panel),  with  the  90%  confidence  bounds  of  the  fit  shown  in  pink.  In  the  right  panel 
of  Fig.  19,  we  have  taken  the  first  and  second  derivative  of  the  PL  signal.  Modeling 
demonstrated  that  the  zero  crossing  (Y  =  0)  of  the  second  derivative  of  the  power- 
dependent  PL  occurs  at  an  applied  field  of  Eo  =  2865  V/cm  when  E0  d12/h  = 
l/-j3T\f£  .  Here,  7)  and  75  are  the  lifetime  of  the  excited  state  and  the  dephasing 
time,  respectively,  and  dn  is  the  dipole  moment  for  the  chosen  2-level  system.  It  is 
significant  to  note  that  the  zero  crossing  of  the  second  derivative  is  within  the 
excitation  range  of  acquired  data  and  not  in  the  extrapolated  portion  of  the  fit. 

The  strong  PL  intensity  observed  in  the  transitions  between  the  4Ig/2  <->•  4F3/2 
manifolds  suggests  that  it  would  make  a  suitable  choice  for  study.  To  determine  the 
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optimal  transition  for  study  as  a  2-level  system,  it  is  first  important  to  note  that  the 
majority  of  the  population  in  the  ground  state  will  begin  in  the  lowest  level  of  the 
4l9/2  manifold.  Also  note  that  the  emission  in  Fig.  1 8  from  the  lowest  level  within 
the  4F3/2  manifold  to  the  lowest  level  within  the  4Ig/2  manifold  is  the  most  intense.  It 
has  been  therefore  identified  that  the  optimal  transition  for  study  as  a  2-level  system 
is  between  the  Ahn  ground  state  to  the  lower  level  in  the  4F3/2  manifold,  which  is  at 
916.2  nm. 


Fig.  19  Left  panel:  Saturation  of  the  916.2-nm  emission  line  using  836-nm  excitation  with 
increased  incident  laser  electric  field  strength.  Right  panel:  First  and  second  derivative  of  the 
best  fit  to  the  saturation  data  in  the  left  panel.  The  value  at  Y  =  0  of  the  second  derivative  is 
related  to  the  dipole  moment  of  the  transition. 


9.3  Characterization  of  Transition  Linewidth 

Demonstration  of  coherence  properties  requires  a  high-quality  sample  with  narrow 
homogeneous  and  inhomogeneous  broadening  (linewidth).  For  this  reason,  the  PL 
linewidth  was  characterized  for  various  growth  conditions.  Figure  20  plots  the  PL 
signal  recorded  for  samples  with  different  doping  levels  and  on  different  substrates: 
E141  at  0.51%  on  sapphire,  E1524  at  0.05%  on  free-standing  GaN  (FS  GaN), 
El 506  at  0.05%  on  sapphire,  and  El 507  at  0.1%  on  sapphire.  The  PL  shows  a 
narrowing  of  linewidth  with  both  a  decrease  in  Nd  doping  and  with  improved 
crystal  quality  on  FS  GaN  substrate  (E1524).  Also,  the  additional  spectra  features 
observed  in  sample  1507,  which  are  attributed  to  Nd  atom  placement  on  minority 
sites  rather  than  the  majority  site,  lead  to  additional  broadening.  While  X-ray 
analysis  has  shown  no  evidence  of  phase  segregation  for  Nd  concentration  below  1 
atomic  %,  it  is  clear  from  the  PL  measurements  that  doping  levels  at  or  below 
0.05%,  grown  on  native  substrates,  are  necessary  to  minimize  spectral  broadening. 
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Linewidth  Comparison 


Wavelength  (nm)  Wavelength  (nm) 

Fig.  20  PL  of  Nd3+:GaN  samples  with  varying  doping  levels.  Right  panel:  Zoomed-in  plot  of 
transition  near  916  nm. 

9.4  Time-Gated  Photoluminescence  (PL)  Setup  and  Experiment 
for  Rabi  Oscillations  and  Photon  Echo  Measurements 

With  the  spectrum  of  Nd3+:GaN  material  for  study  characterized  and  the  transition 
for  the  2-level  system  identified,  the  next  step  involves  time-resolved  experiments 
to  determine  the  lifetime  of  the  state  as  well  as  the  coherence  properties  of  the 
transition  (Rabi  frequency,  photon  echo,  etc.).  For  measurements  where  the 
excitation  laser  and  signal  detection  occur  at  the  same  wavelength,  it  is  difficult  to 
separate  the  2  signals  when  interference  from  scattering  of  the  laser  is  typically 
much  more  intense  then  the  signal  of  interest.  However,  it  is  necessary  to  make 
measurements  in  this  way  in  order  to  observe  transitions  within  and  a  2-level 
system.  We  have  therefore  taken  several  steps  to  help  eliminate  interference  from 
laser  light  scattering.  First,  we  drilled  a  hole  in  the  sample  mount  so  that  the  laser 
may  pass  through  the  sample  and  ensure  that  no  light  is  scattered  off  of  the  cold 
finger.  Second,  polishing  the  backside  of  the  sample  both  reduced  the  overall 
scattering  by  the  sample  and  eliminated  the  scatter  of  randomly  polarized  light. 
Consequently,  the  scattered  laser  light  maintained  its  linear  polarization  and  could 
be  reduced  further  using  a  linear  polarizer  at  the  spectrometer  input.  Finally,  we 
used  the  grating  in  the  spectrometer,  which  diffracts  horizontally  polarized  light 
more  efficiently  than  vertically  polarized  light  by  exciting  the  sample  with 
vertically  polarized  light,  tuning  the  polarizer  to  pass  the  residual  horizontally 
polarized  PL,  which  detected  with  higher  efficiency  because  of  the  polarization 
dependence  of  the  grating.  The  combination  of  these  efforts  reduced  the  detected 
interference  from  laser  scattering  by  a  factor  of  almost  10-4. 

Determining  the  Rabi  frequency  informs  on  the  characteristics  of  the  appropriate 
7i/2  pulse  for  preparation  of  the  Nd3+  ions  into  the  pu  state  and  a  tt  pulse  to  reverse 
the  precession  in  the  wv-plane  for  the  observation  of  a  photon  echo.  For  photon  echo 
measurements,  sample  preparation  occurs  using  a  tt/2  pulse  followed  by  a  n  pulse 
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at  a  variable  time  delay  after  the  n/2  pulse.  Since  this  is  a  coherent  process,  the 
photon  echo  is  emitted  collinearly  with  the  laser  excitation.  Varying  the  time  delay 
between  the  tc/2  and  n  pulses  will  result  in  a  change  in  the  amplitude  of  echo  that 
decays  exponentially  with  the  time  constant  of  the  decoherence  time  72. 

Attempts  to  observe  Rabi  oscillations  were  made  by  preparing  the  excitation  laser 
into  a  square  pulse(s)  of  controllable  and  variable  pulse  width,  and  gating  of  the 
signal  detection  (shown  in  Fig.  21,  left  panel).  The  CW  laser  is  modulated  into 
square  pulses  of  controllable,  varying  widths  between  100  ns  and  100  ps,  using  an 
AOM.  For  the  initial  measurements  we  used  an  electro-optic  modulator  to  chop  the 
laser;  however,  a  2%-5%  laser  radiation  leakage  out  of  the  desired  pulse  window 
clouded  the  results.  Switching  to  modulation  by  the  AOM  improved  the  signal-to- 
noise  ratio  by  more  than  a  factor  of  3.  A  Pockels  cell  was  installed  and  tested  to 
provide  gatable  detection.  This  option  to  the  experimental  system  provides  gated 
detection  with  approximately  4-ns  rise  times  and  windows  between  0  to  1.2  ps.  The 
input  to  the  Pockels  cell  has  been  aligned  with  respect  to  our  fiber  coupling  to 
provide  the  ability  to  switch  between  different  inputs  without  the  need  to  further 
modify  the  system.  A  LabVIEW  program  capable  of  varying  the  width  of  a  single 
pulse,  and  apply  2  pulses  of  varying  width  and  spacing,  was  used. 

Ultimately,  scattering  from  the  sample  still  proved  too  large  to  observe  oscillations 
in  the  PL  emission  at  the  excitation  wavelength.  We  therefore  shifted  our  detection 
wavelength  to  sample  decay  from  the  pi  i  state  to  other  levels  within  the  4  [9/2  manifold 
to  avoid  laser  interference.  This  approach  was  successful  in  removing  effects  from 
scattering;  however,  only  an  approximately  linear  increase  of  the  PL  intensity  and 
not  oscillations  has  been  observed  (Fig.  21,  right  panel).  This  could  be  due  to  a 
number  of  factors,  including  the  relatively  wide  bandwidth  of  the  excitation  laser 
exciting  too  wide  a  bandwidth  of  the  ensemble;  the  sample  temperature  being  too 
high,  leading  to  phonon  broadening;  low  detection  sensitivity/poor  temporal 
resolution;  or  short  coherence  times. 


Fig.  21  Left  Panel:  Timing  diagram  for  the  measurement  of  Rabi  oscillations.  The  gating 
pulse  removed  the  in-pulse  PL  signal  and  laser  scattering.  Middle  Panel:  Sample  time-resolved 
PL  emission  signals  for  a  sampling  of  pulse  widths.  Right  Panel:  Integrated  PL  signal  as  a 
function  of  gating  pulse  width  shows  only  a  linear  increase  in  signal. 
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Realizing  problems  with  our  laser  bandwidth  and  detection  efficiency,  we 
collaborated  with  Prof  Andrei  Faraon  at  the  California  Institute  of  Technology.  The 
Caltech  group  has  been  working  in  the  area  of  quantum  information  of  rare-earth 
doped  solids  for  a  number  of  years  and  is  well  equipped  to  perform  coherent 
spectroscopy  of  Nd-ion-doped  samples.  This  collaborative  work  has  helped  to 
accelerate  the  implementation  of  our  research  plan  and  to  bring  new  experimental 
capabilities  and  understanding  to  ARL. 

During  a  visit  to  Caltech  in  September  2015,  we  characterized  several  in-house- 
grown  samples  at  their  facilities  with  the  intent  of  performing  photon  echo  and  hole¬ 
burning  measurements,  and  to  try  quantum  memory  experiments  based  on  AFC 
protocol.  Using  their  frequency-locked,  narrow  linewidth  laser,  we  were  able  to 
observe  an  approximately  80-pm  emission  linewidth  to  the  PL  signal  near  916  nm 
when  a  higher  manifold  was  excited,  demonstrating  a  good  homogeneous 
linewidth.  The  PL  decay  measurement  yielded  a  result  of  T\  =  150  ps  (Fig.  22),  but 
we  did  not  observe  a  photon  echo.  There  are  2  reasons  why  we  were  not  able  to 
observe  a  photon  echo:  1)  the  signal  is  very  weak  or  2)  Ti  is  very  short.  Also,  the 
exponential  decay  of  the  photon-echo  amplitude  at  a  characteristic  time  of  Ti 
implies  that  a  short  T2  will  lead  to  a  weak  signal.  The  weak  signal  for  our  samples 
is  likely  due  to  the  thin-film  growth  (approximately  hundreds  of  nanometers), 
which  provides  orders  of  magnitude  less  material  (and  interaction  length)  than  the 
approximately  0.5-mm  samples  the  Caltech  group  typically  studies.  To  address  this 
problem,  a  waveguide  is  necessary  to  increase  the  interaction  length  of  the 
interrogation  laser.  Modeling  has  demonstrated  that  it  is  not  possible  to  grow  a 
coherently  strained,  quality  Nd-doped  GaN  film  of  sufficient  thickness  to  confine 
the  laser  mode  on  a  layer  of  material  with  a  satisfactory  index  mismatch.  Two 
alternative  approaches  were  identified  to  attain  mode  confinement.  First,  the  GaN 
substrate  could  be  removed  using  a  process  such  as  photoelectrochemical  (PEC) 
etching,  which  would  target  a  sacrificial  layer  of  narrower  bandgap  grown  between 
the  Nd-doped  GaN  layer  and  the  GaN  substrate.  Then,  waveguide  structures  could 
be  etched  out  of  the  remaining  Nd-doped  GaN.  Second,  an  (Al)GaAs  layer  could 
be  grown  and  etched  into  a  waveguide  on  top  of  a  high-quality  Nd3+:GaN  film  that 
would  take  advantage  of  the  coupling/absorption  of  the  evanescent  wave  into  the 
Nd3+:GaN.  These  approaches  would  require  significant  effort  and  development 
outside  of  this  program. 
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Fig.  22  Measurement  of  T\  =  150  ns  in  Nd3+:GaN 

9.5  Characterization  of  Field-Shifting  of  Transition  Wavelength 

In  parallel  with  experiments  on  the  coherence  properties  of  Nd3+:GaN,  we  also 
investigated  the  effect  of  varying  the  electric  field  on  the  Nd  PL  spectrum  to  provide 
a  tunable  memory.  To  vary  the  applied  field,  we  designed  and  grew  a  series  of  Nd- 
doped  GaN  p-i-n  structures,  strain-balanced  superlattice  (SBSL)  structures,  and 
periodic  unbalanced  (PUB)  structures.  In  the  SBSL  and  PUB  structures,  the 
spontaneous  polarization  discontinuity  between  AlGaN  and  GaN  layers  enable  a 
large  range  of  electric  fields  to  be  generated  along  the  c-direction.  Modeling  of  the 
built-in  electric  fields  was  performed  using  a  self-consistent  drift-diffusion  model. 
In  the  SBSL  design,  higher  fields  were  targeted  using  periodic  structures  grown 
with  alternating  layers  of  GaN  and  AlxGai-xN  with  equal  thicknesses,  strain- 
balanced  to  a  template  layer  of  Al^Gai-^N.  The  PUB  structures  were  necessary 
to  target  fields  below  approximately  500  kV/cm  because  of  the  large  thicknesses 
required  in  the  SBSLs  for  low-field  intensities,  which  would  likely  lead  to  misfit 
defect  formation  and  increase  the  inhomogeneous  broadening  of  the  Nd3+  ion 
transitions.  This  was  accomplished  using  thin,  periodic,  low-Al  mole  fraction 
AlGaN-strained  layers  grown  pseudomorphically  to  a  free-standing  GaN  substrate 
to  induce  a  field  in  the  alternating  Nd-doped  GaN  layers. 

Figure  23  summarizes  the  modeled  electric  field  ranges  possible  for  different 
AlGaN/GaN  SBSL  and  PUB  designs.  We  have  completed  the  layer  structure  designs 
for  both  the  SBSL  and  the  PUB  structures,  and  we  have  grown  structures  with 
calculated  electric  field  strengths  ranging  from  274  kV/cm  to  4.5  MV/cm  by  growing 
2  series:  SBSL  structures  using  Alo.8Gao.2N  barriers  on  an  Alo.4Gao.6N  template  layer 
(orange  line)  and  the  PUB  structures  with  Alo.2Gao.8N  barriers  on  a  FS  GaN  substrate 
(green  line).  SIMS  measurements  were  performed  to  validate  that  Nd  incorporation 
occurred  in  the  desired  layer.  Figure  24  plots  the  SIMS  measurement  on  the  60-nm- 
thick  well/barrier  SBSL  sample  that  demonstrates  that  Nd  incorporation  is  present  in 
the  GaN  layers  (within  the  resolution  of  the  measurement). 
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Fig.  23  Electric  field  vs.  GaN  well/  AlxGa(i-X)N  barrier  thickness  for  strain-balanced 
superlattice  (SBSL)  structures  with  GaN  wells/  Al0.sGa0.2N  barriers  (orange  line)  and  periodic 
unbalanced  (PUB)  structures  with  GaN  wells  and  Al0.2Ga0.sN  barriers  (green  line) 


E1518 


Fig.  24  SIMS  measurement  on  the  60-nm-thick  well/barrier  SBSL  sample 

Figure  25  plots  the  PL  spectra  of  the  series  of  high-field  SBSL  samples  (with 
designed  electric  fields  of  640  kV/cm,  1.0  MV/cm,  2.3  MV/cm,  and  4.5  MV/cm) 
and  an  approximately  zero-field  600-nm  bulk  Nd-doped  GaN  sample  and  a  high- 
field  (-4.5  MV/cm).  PL  spectra  were  acquired  using  a  tunable  pulsed  Ti:Sapphire 
laser  with  wide  spectral  bandwidth  tuned  to  excite  transitions  from  the  4Ig/2  ground 
state  to  the  4F5/2  manifold  near  836  nm  at  10  K  (to  ensure  that  most  of  the  population 
in  the  ground  state  would  begin  in  the  lowest  level  of  the  4I9/2  manifold).  PL  from 
the  lower-energy  4F3/2  manifold  to  the  Aha  manifold  was  detected  using  a  Fourier- 
transform  infrared  spectrometer  (FTIR)  to  separate  transitions  and  capture  the 
entire  spectrum  at  once.  This  broadband  approach  was  employed  to  mitigate  the 
difficulty  that  results  from  excitation  and  emission  linewidths  shifting  spectrally.  A 
nearly  identical  shift  of  greater  than  1 . 1  nm  is  observed  in  the  PL  spectra  for  each 
sample,  demonstrating  the  ability  to  tune  the  Nd  transition  with  field.  However,  the 
magnitude  of  the  shift  observed  (>300  GHz)  is  more  than  an  order  of  magnitude 
larger  than  expected  for  a  Stark  shift  at  these  built-in  electric  fields. 
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Fig.  25  Nd3+:GaN  PL  spectra  of  an  approximately  zero-field  bulk  sample  and  high-field 
SBSL  samples  (expected  electric  field  strength  given  in  legend) 

We  therefore  began  a  collaboration  with  Profs  David  Smith  and  Martha  McCartney 
at  Arizona  State  University  to  perform  electron  holography  on  our  samples  to 
provide  a  direct  measurement  of  electric  field.  Results  from  that  study  were  both 
informative  and  perplexing.  Transmission  electron  microscopy  (TEM)  images  are 
shown  in  Fig.  26  for  sample  E1521;  the  40-nm-thick  well/barrier  SBSL  sample  is 
expected  to  possess  a  built-in  electric  field  of  approximately  1  MV/cm.  It  is 
apparent  that  the  layered  structure  is  compromised  with  irregular  interfaces  and 
facets  that  are  likely  due  to  propogation  of  defects  from  the  substrate  layer  (see 
upward  “lines”  in  the  base  of  Fig.  26,  left  panel,  that  seem  to  lead  to  sharp  structures 
farther  up  in  the  sample).  The  propogation  of  these  defects  likely  expose  different 
crystal  planes  (rather  than  the  desired  c-plane),  which  posses  differing  growth  rates 
and  lead  to  the  observed  structuring  of  the  periodic  sample. 


Fig.  26  TEM  images  of  sample  1521,  the  40-nm-thick  well/barrier  SBSL  sample,  with  a  100- 
nm  scale  bar  (left  panel)  and  a  20-nm  scale  bar  (right  panel). 

The  TEM  images  of  sample  1521  were  processed  using  an  electron  holography 
approach  described  in  McCartney  and  Smith.14  Briefly,  the  technique  images  the 
phase  shift  of  an  electron  beam  that  passes  through  a  sample  with  one  that  passes 
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through  vacuum.  A  sample  phase  image  for  sample  E1521  is  shown  in  Fig.  26  (left 
panel).  The  brighter  areas  are  GaN,  where  the  built-in  electric  field,  E,  is  indicated, 
and  the  darker  areas  are  AlGaN.  The  area  within  the  outlined  area  in  Fig.  26  (left 
panel)  is  plotted  as  a  line  profile  of  phase  (in  radians)  in  Fig.  26  (right  panel).  The 
formula  for  converting  the  phase  image  to  electrostatic  potential  is 

A(p(x,  y)  =  f*  CEV (x,  y,  z)dz~  CEV (x,  y)t(x,  y)  ,  (30) 

where  A cp  is  the  phase  difference  (from  vacuum,  for  example,  or  from  one  side  of 
the  layer  to  the  other),  Ce  is  a  constant  that  depends  on  accelerating  energy,  and  t 
is  the  thickness.  Equation  30  assumes  that  the  potential  is  constant  throughout  the 
thickness,  which  ignores  the  surface  effects  in  the  thin  film,  but  since  the  sample  is 
approximately  tens  to  hundreds  of  nanometers  thick,  “bulk”  effects  can  be  assumed 
to  dominate.  Another  important  approximation  assumes  that  there  are  no  other 
phase-shifting  mechanisms.  Estimates  of  the  built-in  electric  field  throughout  the 
sample  (beyond  the  shown  images  area  in  Fig.  27)  yield  electric  fields  in  the  AlGaN 
layers  on  the  order  of  0.10  MV/cm,  and  the  field  in  the  GaN  layers  to  be  in  the 
0.10-0.25  MV/cm  range.  Each  estimate  is  considerably  lower  than  what  is  modeled 
for  this  structure,  where  an  approximately  1 -MV/cm  electric  field  magnitude  is 
expected  in  the  GaN  layer. 


Substrate 


Fig.  27  Electron  holography  phase  image  (left  panel)  with  area  for  line  profile  (right  panel) 
outlined.  Brighter  areas  are  GaN,  and  darker  areas  are  AlGaN. 

The  lower-than-expected  estimate  further  complicates  the  larger-than-expected 
spectral  shift  of  the  PL  wavelength  for  the  identified  transition.  The  clear 
implication  is  that  the  spectral  shift  is  due  to  some  other  effect  than  the  built-in 
electric  field  due  to  the  heterostructure.  An  initial  theory  that  the  overall  crystal 
field — due  to  the  on-average  crystal  structure  of  Alo.4Gao.6N,  coherently  strained  to 
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a  base  layer  of  Alo.4Gao.6N — is  dominating  the  transition  wavelength  has  been 
eliminated  since  the  transition  in  the  binary  AIN  shifts  to  shorter,  not  longer, 
wavelengths  than  the  Nd3+:GaN  transition.  A  remaining  theory  suggests  the 
mechanism  that  shifts  the  PL  wavelength  is  due  to  the  complex  strain  state  of  the 
highly  defective  sample,  which  results  in  space  charge  at  the  defect  locations. 
Further  investigation  of  this  phenomena  will  focus  on  samples  grown  on  high- 
quality  substrates,  such  as  the  set  of  low-field  PUB  samples  and  p-i-n  structures 
that  can  be  biased  to  determine  the  functional  dependence  of  the  shifted  Nd 
transition  on  applied  field. 

10.  Conclusion 


In  this  DRI  program,  we  investigated  the  potential  of  a  cryogenically-cooled  REIC, 
Nd3+:GaN,  for  use  as  tunable  solid  state  quantum  memory.  Theoretically,  we 
developed  a  quantum  description  of  2-  and  3 -level  systems  and  of  the  2-  and  3- 
pulse  photon  echo  effect,  and  estimated  the  effect  of  an  applied  electric  field  for 
control  of  collective  coherences  using  the  linear  Stark  effect.  Experimentally,  PL 
measurements  were  used  to  study  the  energy  level  structure  of  the  Nd  ions 
embedded  into  the  crystal  lattice  of  GaN.  The  power  dependence  of  the  PL  emission 
intensity  was  used  to  determine  a  relation  between  key  parameters,  while  the  decay 
rate  of  the  PL  emission  yielded  an  excited  state  lifetime  of  7W50  ps.  SBSL 
structures  with  high  built-in  fields  demonstrated  the  ability  to  shift  the  Nd  transition 
wavelength,  however  results  were  dominated  by  sample  defects.  Additional  work 
is  required  to  determine  the  dependence  of  the  Stark  shift  on  electric  field  and  to 
measure  the  coherence  times  in  Nd3+:GaN  waveguide  structures. 
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List  of  Symbols,  Acronyms,  and  Abbreviations 


3-D 

AFC 

AIN 

AOM 

ARL 

at% 

CRIB 

CW 

DC 

DRI 

FS  GaN 

FTIR 

GaN 

He 

MBE 

Nd 

OPA 

PA-MBE 

PEC 

PL 

PLE 

PUB 

REIC 

RF 

RS 


3-dimensional 
atomic  frequency  comb 
aluminum  nitride 
acousto-optic  modulator 
US  Army  Research  Laboratory 
atomic-percent 

controlled  reversible  inhomogeneous  broadening 
continuous  wave 
direct  current 

Director’s  Research  Initiative 
free-standing  gallium  nitride 
Fourier-transform  infrared  spectrometer 
gallium  nitride 
helium 

molecular  beam  epitaxy 
neodymium 

optical  parametric  amplifier 
plasma-assisted  molecular  beam  epitaxy 
photoelectrochemical 
photoluminescence 
photoluminescence  excitation 
periodic  unbalanced 
rare-earth-ion-doped  crystal 
radio  frequency 
Russell-  S  aunders 
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RWA 

rotating-wave  approximation 

SBSL 

strain-balanced  superlattice 

SIMS 

secondary  ion  mass  spectrometry 

TEM 

transmission  electron  microscopy 

XRD 

X-ray  diffraction 
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